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We calculate the power density required to excite various states of 2 and which may 
be expected to radiate light observable from the ground. These calculations are made by 
first calculating the electron energy distribution functions in the presence of the KF electric 
fields, and from this the rate of excitation of the radiating states. We discuss the feasibility 
of experiments to generate sufficient intensity to cause appreciable modification in the night 
airglow intensity. The possibility of performing experiments using the two body dissociative 
attachment of electrons is also considered, especially in regard to the performance of basic 
experiments in atomic physics. 



1. Introduction 

In the discussion of nonlinear processes in the 
ionosphere, the excitation of optical radiation by a 
radio beam is an extreme case. Yet, if transmitters 
sufficiently powerful to excite optical radiation in the 
ionosphere existed, a great deal of quantitative 
information concerning the constitution of the iono- 
sphere could be collected on a routine basis. It is, 
therefore, of some importance to investigate the 
conditions under which detectable optical radiation 
might be excited. 

This problem has been considered by others, 
Bailey [1959], Clavier [1961]; but not in the detail 
which we will apply. Authors in the past have esti- 
mated the intensity to be expected in the atmosphere 
based on laboratory data obtained in a discharge. 
Such a procedure is liable to considerable error, since 
conditions in a discharge differ greatly from those 
extant in the atmosphere. In particular, in dis- 
charges there exist many metastable states which 
change the power required to heat electrons and 
recombination processes which yield a great deal of 
the luminosity. In addition, the electron density in 
discharges is often such that the electron energy dis- 
tribution is dominated by electron interactions, a 
condition which does not apply in the lower iono- 
sphere, [Megill and Carleton, 1964]. We shall, in 
this paper, concern ourselves principally with the 
direct "heating" of the electrons already in existence 
in the ionosphere by means of electric fields and 
with the excitation of electronic states in the neutral 
gas by collisions with these electrons. The problem 
of "breaking down" the atmosphere will be discussed 
only briefly. 



2. Excitation of Optical Radiation 

If we know the distribution i(u) of electrons with 
energy, u, under the influence of electric and magnetic 
fields then we can calculate the total rate of excitation 
vy of a particular state from 
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where N t is the number density of the gas, N e is the 
number density of electrons, tr# is the cross section 
for excitation to state j, v is the electron velocity, 
and we have normalized the distribution function 
so that 
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Carleton and Megill [1962] have described the 
calculation of the electron energy distribution func- 
tion for altitudes of interest here. It is shown that 
the distribution function may be determined by 
solving the following differential equation. 
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where 



k is the Boltzmann constant. 

T is the temperature of the neutral particles. 

v is the elastic collision frequency 2 N t at v where 
Ni is the density of gas i and a is the elastic scattering 
cross section for momentum transfer. 
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6 is the average fractional energy lost per elastic driving function e by which the electrons gain energ 
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collision -jj where m is the electron mass and M is 

the mass of the gas particle. 

7 is — =-( — ) e being the electronic charge and E 
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the mis electric field, 

Q{u) is 2 A^- 3/2 (m/2)" 1/2 x 
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where o- and JV* are as in (1) and u {j is the energy 
required to excite the state u is . The summation is 
overall inelastic states j for each gas i. 

a is a scaling parameter on the electric field 
expressed as 

, 2 +co 2 +col 
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where co is the angular frequency of the exciting 
field E and co B is the gyrofrequency of the electron in 

eB 

the earth's magnetic field B where w B = : — 

m 

The solution of (3) is carried out numerically in 

order to obtain the desired distribution function, 

using experimental cross sections where available, and 

theoretically estimated ones elsewhere. A list of the 

cross sections used for these and other calculations 

has been given by Megill and Carleton [1964]. 

3. Scaling Parameters 

While one may calculate the expected excitation 
for any given power density, excitation frequency, 
and collision frequency, the number of such calcula- 
tions required to cover all possible cases is very large. 
We search, therefore, for scaling parameters in the 
theory which will allow more general application of 
calculated values. 

If we investigate (3) in detail, we find that (after 
dividing through by v(u) for convenience) the left- 
hand side of the equation constitutes a kind of 
"driving function" which describes the manner in 
which the electrons gain energy while the right-hand 
side describes the mechanisms which will cause the 
electrons to lose energy. We further see that on the 
right hand side of (3) only the parameters of the gas 
enter. The only reference to the mechanism which 
increases the electron energy then, is inside the 
brackets on the left hand side of the equation 
(kTd+ya). Recall that we have divided the equa- 
tion by v(u). 

The term kT8 presents a base energy (due to the 
neutral gas temperature) which is the electron 
energy for zero electric field. We have, then, as the 
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Now let us examine expression (4) for limiting cases. 
In the ionosphere we have a rapidly decreasing 
number density with height so that the collision 
frequency v = Nav is also a rapidly varying function 
of height. Consider a beam of radiation extending 
upward through the ionosphere. It will, except for 
a small transition region, be in a situation in which 
v^>^>oo-\-q)b or v<^<^o) B . An exception is the con- 
dition co = co s about which more will be said later. 
We will, accordingly, consider these two cases in 
sequence. 



When v<^ co and if co s <^co 
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which indicates that the mechanism by which the 
electrons gain energy is independent of the pressure 
of the gas and is only dependent upon its character 
as described by the right-hand side of (3). 

On the other had, if v^> co ± w B then (4) reduces to 
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which, since v oc N and N is proportional to the 
pressure, yields the scaling factor E/P which is 
commonly used in gaseous electronics. 

When co = co B , gyroresonance exists. For pur- 
poses of calculating distribution functions, this 
may be considered to be the same as if there were a 
d-c field. This is strictly true only if the radiation 
is circularly polarized in the extraordinary mode. 
For most purposes of interest in ionospheric prob- 
lems, however, such an approximation is adequate. 

In general, the electric field of importance in 
these calculations will be that present in the iono- 
sphere at the height for which the calculations are 
being made. This field is, of course, dependent 
upon the power of the transmitter, the gain of the 
antenna and the absorption of the radiation below 
the point of interest. The calculation of the 
absorption of the radiation is in itself a difficult 
problem which is not discussed here. 

We shall define two parameters then, in terms 
of which we shall present the data. These we will 
call the HF or high frequency and the LF or low 
frequency reduced power (P r ) H F and (P r ) L F where 



and 
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If E is in volts per cm and n is in ohms then these 
quantities have units of watts cm" -2 sec 2 and watts 
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cm~ 2 cm 6 or watts cm 4 if N is in units of number per 
cubic centimeter. These units, while they may 
appear unwieldy, are convenient for application 
to geophysical problems. 

The power required from an antenna in order to 
generate a specific effect may in turn be calculated 
from 
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where h is the height at which the excitation rate 
is required, N is the number density at that height 
and G is the gain of the antenna. T is the amount 
of power transmitted to height h. For the purposes 
of our further calculations T is considered to be one. 



4. Observable Radiations 

Before calculating excitation rates of various 
states, let us first consider what radiations are 
likely to be visible from ground level. The first 
condition, of course, is that the experiment be done 
at night. Because we are dealing here with excita- 
tion by low energy electrons which already exist in 
the ionosphere as opposed to those created in a 
discharge we need consider radiation from states 
which can be excited from ground state atoms or 
molecules. The density of metastable states which 
can be excited to higher levels and will then radiate 
are negligible. Ground state rotational or vibra- 
tional states of the major molecular constituents 
will not radiate, and so, energy put into these states 
will not result in visible radiation. We discuss here 
only the excitation of electronic Levels in the atoms 
and molecules which may possibly reradiate. In 
calculation of the distribution function, however-, all 
inelastic states have been considered even though 
they are not discussed explicitly in this paper. In 
nitrogen, such states require so much energy to be 
excited (greater than 6.2 eV) that excitation is 
unlikely except in extremely powerful systems. In 
2 the three most easily excited electronic states are 
the a l L g (0.98 eV), the &%+ (1.65 eV), and the 
A 3 2 M + (4.47 eV). Of these, the first two have 
extremely small cross sections for direct excitation 
by electron collisions [Schulz and McDowell, 1962], * 
so that we need to discuss only the excitation of the 
A 3 2 M + state which radiates in the Herzberg bands. 
These bands are seen in the night airglow and are 
rendered difficult to observe because the energy is 
distributed into many lines in the band. However, 
fairly small amounts of excitation should be 
detectable with sophisticated techniques. 



This leaves the X D (1.96 eV) and l S (4.17 eV) 
states of atomic oxygen. The former would be 
undetectable because it is strongly deactivated by 
collisions with 2 at altitudes below about 140 km. 
These calculations do not apply in any practical 
situation for altitudes above 140 km because the 
dominant energy loss term there would be that due to 
diffusion out of the exciting beam, an effect which is 
not included in the calculation of the distribution 
function. 0( X S) is the strongest line of the normal 
night airglow and variations of intensity should be 
readily observed if excitation is accomplished. 
This leaves then, as the two most readily observed 
emissions 2 (A 3 2 M + ) and 0( X S) both of which require 
more than 4 eV for excitation. The generation of 
electrons of this energy is particularly difficult in air 
because there is a very large cross section for energy 
loss to the vibrational states of N 2 for all electrons 
between about 1.7 and 4 eV [Schulz, 1959]. It is of 
interest to know, however, just what power would be 
required in order to obtain detectable intensities. 
In figures la and It), we plot specific rates of excita- 
tion for the above two states, as well as two body 
attachment rates and ionization rates, as a function 
of the two previously defined reduced powers. In 
order to obtain the excitation rate at a particular 
power and height, one needs to multiply the value 
given in the figure by electron density and constituent 
density. The cross section for the excitation of the 
2 (A 3 2 W +) is estimated from the data of Schulz 
[1962], that for the excitation of 0( ! S) from the 
theoretical calculations of Seaton [1953]. As an 
example let us suppose that we wish to know the 
power required to excite 100 rayleighs 2 of 5577 A 
radiation (from ()( 1 S) in a slab of the ionosphere 
between 90 and 100 km). Then this means that we 
must excite 100 atoms/cm 3 /sec, which will later 
radiate. If we assume that there are on the average 
500 electrons/cm 3 and 2X10 12 oxygen atoms/cm 3 
in t he slab, then we would expect to see approximately 
100 rayleighs for a specific excitation rate of 0( X S) 
equal to \0~ u which requires from figure la a reduced 
HF power of 3.3X10" 21 W sec This results (using 
5a) in a required power flux of 8X10 -5 W/cm 2 , if 
we are considering a transmitter at 50 Mc/s. Using 
(6a) it is found that 9X10 10 W are required for an 
isotropic radiator. If we have an antenna gain of 
10 4 , we then need 9 MW transmitted from the antenna. 
The very large power required for minimum detec- 
tion conditions indicates that experiments resulting 
in the excitation of optical radiation may not be 
economically feasible unless the installation can be 
justified for other applications. 

An additional consideration should be made at this 
point. The calculations made here yield average 
excitation rates while the power is turned on, and 
after a steady state has been established. The time 
required to establish equilibrium is approximately 
equal to the product of the average fractional energy 
loss per collision (averaged over all types of collisions) 



1 These states are, of course, excited to some extent; but with the method of 
excitation considered here, we find it will never, at any power less than that 
which is required for breakdown, be possible to excite a detectable amount of 
radiation from them. 



2 The rayleigh is defined as Iff 5 emissions per cm 2 column per second through 
the radiating medium. 
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Figure la. This figure shows (1) the rate of excitation ofO ( J D) , 
(2) the rate of attachment to 2 , (8) the rate of excitation of 
2 (A 3 2 U + ), and (4) the rate of ionization as a function of 
reduced power with high frequency excitation. 

To obtain an absolute rate the ordinate should be multiplied by the constituent 
density and the electron density. The power flux in watts cm- 2 may be had by 
multiplying the abscissa by (2tt/) 2 where/ is the frequency of excitation. These 
curves should be used if v <co±w b . 




Figure lb. This figure shows the same quantities as figure 
la except that they are plotted versus the reduced power for low 
frequency excitation. 

In this case the abscissa should be multiplied by the square of the total number 
density in particles per cm3 to yield watts cra-l These curves should be used 
when j>>o±w& and for the gyroresonance condition. 



and the collision frequency. For air the average 
energy loss is not constant with energy, but for 
rough calculations may be assumed to be of the order 
of 10~ 3 . For the convenience of the reader and 
those who might wish to extend the calculations to 
include absorption we plot in figures 2a and 2b the 
average collision frequency, normalized to one par- 
ticle/cm -3 and the average energy against the two 
previously defined reduced powers. The relaxation 
time will vary from the order of a microsecond at 60 
km height to the order of a millisecond at 100 km 
height. Transmitters must, therefore, be operated 
with pulse lengths long compared to these times for 
steady state to be achieved. For the higher heights 
this puts severe requirements on the transmitter 
because of the long pulses required. 

One other parameter which must be considered 
for the case of 0( 1 S) is that the state has a mean life 
of the order of 1 sec. For this reason the effect of 




(Pr) hf 

Figure 2a. The abscissa is as defined in figure la while the 
ordinates are the average collision frequency {curve 2) per 
neutral particle in sec~ l and the average energy in eV (curve 1). 
The averages are over the calculated distribution functions. 




(Pr), 



Figure 2b. The abscissa is as defined in figure lb while the 
ordinates are as defined in figure 2a. 
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pulsed operation will be to reduce the observed 
change in intensity by the duty cycle. It is to be 
noted that this does not mean that the average 
power flux is the important quantity. The ex- 
tremely nonlinear nature of the intensity versus 
power curve indicates that the transmitter should be 
operated at as high a peak power as possible and then 
for as high a duty cycle as is compatible with this 
condition. 

5. Modification of the Ionosphere by 
Attachment 

Another possible experiment of interest is the modi- 
fication of the atmosphere by attaching the electrons 
to 2 molecules. The problem has been discussed in 
some detail by Molmud [1964] for attachment by 
three body attachment, 



2 +e+M^0 2 -+M. 



(7) 



The rate of this process is probably the easiest one 
to change with transmitted RF power because there 
is a maximum in the rate at about 0.1 eV average 
electron energy [Chanin, Phelps, and Biondi, 1959]. 

There is a second process which occurs at higher 
energies. This is dissociative attachment of 2 , by 
the following: reaction, 



O 2 +e-^Cr+0. 



(8) 



This process has a threshold at 4 eV with a maximum 
at 6.3 eV [Schulz, 1962]. We may calculate this 
rate in the same manner that one calculates attach- 
ment rates. These rates are given in figures la and 
lb. In the example previously discussed, if we were 
exciting 100 rayleighs of the 5577 A line we would 
expect electrons to attach at a rate of 3 percent per 
second. The electrons would effectively disappear 
at this rate until they are re-released either by colli- 
sional detachment or by the reverse of process (8). 
These rates are not well known and it might be that 
an ionospheric experiment would be one of the best 
with which to determine their rates. 

6. Breakdown in the Ionosphere 

The general technique of calculation used here 
can, in principle, be used to compute the power 
required to "breakdown" the ionosphere. In the 
lower ionosphere one may, to a large extent, ignore 
diffusion losses and may therefore compute the rate 
of attachment as a function of power and the rate 
of ionization as a function of power. As the ioniza- 
tion rate becomes the larger the ionization density 
will grow exponentially until some limiting process 
sets in. The details of an approximate treatment of 
some of these processes is given by Lombardini 
[1964]. In addition he gives comparisons to experi- 
mental data which indicate the order of the accuracy 
of the calculations. Specific rates of ionization as a 
function of {P t )hf and (P t )lf are shown in figures 
la and lb. 



7. F-Region Emissions 

In the investigation of possible modification of 
airglow intensities another possibility of interest pre- 
sents itself in terms of the /-region emission from 
0( 1 D) at 6300 A. This radiation is presumably due 
to the two step recombination of + by means of 



and 



0++N 2 -^NO++N 



NO+^-c-^N+OCD). 



(9) 



(10) 



The rate of (9) is the slower and controjs the emission 
rate so that the intensity of 6300 A radiation is 
proportional to the product of the + and N 2 
densities. 

Farley [1963] has shown that by using power at 
the plasma frequency it is possible with reasonable 
power to heat (—250 kW) the electrons enough 
that the enhanced diffusion decreases the electron 
and ion density by about 10 percent or more. This 
should, in turn, decrease the intensity of 6300 A ra- 
diation by 10 percent and serve as an experimental 
check on the recombination process in the F region. 
It should be noted that the mean life of the 0( 2 D) 
state is 100 sec, so that a transmitter which will 
perform the described experiment must be operated 
CW. 

8. Conclusions 

The requirements on large power transmitters 
which can be expected to yield observable amounts 
of optical radiation have been expressed quantita- 
tively for the two most likely radiations. We find 
that the size of the required installations is so large 
that such experiments may not be economically 
feasible at the present time. 



The author recognizes the very considerable con- 
tributions of Mrs. J. D. Droppleman in the perform- 
ance of the calculations and the preparation of the 
figures. 
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